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For many-particle systems defined on lattices we investi- 
gate the global structure of effective Hamiltonians and ob- 
servables obtained by means of a suitable basis transforma- 
tion. We study transformations which lead to effective Hamil- 
tonians conserving the number of excitations. The same 
transformation must be used to obtain effective observables. 
The analysis of the structure shows that effective operators 
give rise to a simple and intuitive perspective on the ini- 
tial problem. The systematic calculation of n-particle irre- 
ducible quantities becomes possible constituting a significant 
progress. Details how to implement the approach perturba- 
tively for a large class of systems are presented. 

PACS indices: 75.10.-b, 75.40. Gb, 03.65.-w, 02.30. Mv 



I. INTRODUCTION 

Effective models are at the very center of theoretical 
physics since they allow to focus on the essential physics 
of a problem without being distracted by unnecessary 
complexity. Hence it is very important to dispose of 
systematic means to derive effective models. Here we 
will present the mathematical structure of a certain kind 
of effective models, namely effective models where the 
elementary excitations above the ground state can be 
viewed as particles above a complex vacuum. This type of 
view is very common in low-temperature physics. Many 
experiments can be understood on the basis of this pic- 
ture. 

In this paper, we will elucidate the global structure of 
the Hamiltonian and of the observables if the model is 
transformed to a model which conserves the number of 
particles. Such a mapping is often possible and renders 
the subsequent calculation of physical quantities much 
easier. The determination of the effective Hamiltonian 
is facilitated by the decomposition into n-particle irre- 
ducible parts. We set up such a classification at zero 
temperature for strong- coupling situations, i.e. no weak- 
coupling limit is needed and no non-interacting fcrmions 
or bosons are required. Generically, we deal with hard- 
core bosons. 

The necessity for the decomposition into n-particle ir- 
reducible parts has arisen in perturbative calculations of 
the effective Hamiltonians because only the n-particle ir- 
reducible interactions are independent of the system size. 
The second main point of this article is the perturbative 
computation of effective Hamiltonians and observables. 
Such computations are a standard technique for ground 



state energies (0-particle terms) and dispersion relations 
(1-particle terms), see Ref. [1] and references therein. 
But the possibility to compute multi-particle contribu- 
tions has only recently been realized [2-4] and continues 
to be exploited intensively. The key ingredient is to de- 
fine a similarity transformation on the operator level (see 
below) . 

A promising alternative route, which we can only 
sketch in this article, consists in the non-perturbative, 
renormalizing realization of the transformation of the ini- 
tial model to the effective model which conserved the 
number of particles. Examples of this approach are real- 
ized in fermionic models [5-7]. 



A. Starting point 

We consider models which are defined on a lattice T. 
At each site of the lattice the system can be in a number d 
of states spanning the local Hilbert space. Let us assume 
that d is finite. The dynamics of the system is governed 
by a Hamiltonian H acting in the tensor-product space 
of the local Hilbert spaces. For simplicity we do not con- 
sider antisymmetric, fermionic situations although this 
is also possible. So we are focusing on physical systems 
which can be described in terms of hard-core bosons. 

The Hamiltonian H is assumed to be of finite range. 
This means that it is composed of local operators h u act- 
ing on a finite number of sites in the vicinity of the site 



H 



2> 



We further assume, that H can be split as 
H[x) = U + xV , 



(1) 



(2) 



so that the spectrum of U is simple (see below) and that 
the system does not undergo a phase transition from x = 
to the range of values we are finally interested in. These 
requirements do not necessarily imply that x has to be 
small. But it is helpful if this is the case. 

The ground state of U and its lowest lying eigen-states 
shall be known. The latter will be viewed as elementary 
excitations from which the whole spectrum can be built. 
We assume that we can view the elementary excitations 
above the ground state as (quasi-)particles above the vac- 
uum. For simplicity, we will drop the prefix 'quasi-'; it 
is understood that 'particle' is a synonym for elementary 
excitation. 
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We assume that the physical picture sketched for 
H(x = 0) = U is linked continuously to the range 
< x < x c where x c is the critical value at which a 
phase transition occurs. At the critical value x c the pic- 
ture breaks down and cannot be used beyond x = x c . 
Generically, a mode of H (x) will become soft at x c . 

Furthermore, the particles for x = shall be local in 
the sense that we can assign a site to each of them. Let 
Q be the operator that counts the number of particles. 

As a concrete example, the reader may think of an an- 
tifcrromagnetic Hcisenbcrg model made up from strongly 
coupled (coupling J) pairs of spins ('dimers') which 
are weakly coupled (coupling xJ) among themselves, 
e.g. [3,8]. At x = 0, the ground state is the product state 
with singlets on all dimers; the elementary excitations 
are local triplets. The number of these local triplets, i.e. 
the number of dimers which are not in the singlet state, 
shall be given by the operator Q. 

A considerable simplification of the problem can be 
achieved by mapping the initial problem H (x) to an ef- 
fective Hamiltonian H b q{x) in which the number of ele- 
mentary excitations does not change. That is the num- 
ber of particles should be a conserved quantity. Then the 
computation of many physical quantities is significantly 
simplified. 

In this article, we advocate to use a continuous uni- 
tary transformation (CUT) [9-11,8] in order to achieve a 
systematically controlled mapping of the kind described 
above which leads to 

[H eS ,Q}=0, (3) 

i.e. iJ ff conserves the number of particles. Such an ap- 
proach has three major advantages: 

1. Conceptual clarity 

Using a unitary transformation guarantees that no 
information of the orginial model is lost. In par- 
ticular, it is clear that the same transformation 
[9,12-14] can be applied to obtain the effective ob- 
servables O e ff from the original observables O. 

2. Technical simplicity 

To implement the unitary transformation in a con- 
tinuous fashion only the computation of commuta- 
tors is required since the mapping is split into in- 
finitesimal steps leading to a differential equation 
[9] 

d e H=[ V (£),H(£)} (4) 

where I G [0, oo] is an auxiliary parameter 
parametrizing the continuous transformation with 
starting point I = and end point I = oo. 

3. Good controllability 

By an appropriate choice of the infinitesimal gen- 
erator 7] of the transformation it can be designed 
such that is preserves block-band diagonality [8,11]. 



Moreover, it is renormalizing in the sense that ma- 
trix elements between energetically very different 
states are transformed more rapidly than those be- 
tween energetically adjacent states [5,6,11]. 

We like to point out, however, that the general struc- 
ture of operators does not depend on the details of the 
method by which the effective particle-conserving model 
is obtained. Also other methods than CUTs are conceiv- 
able. 

In the present paper we will focus on perturbative re- 
alizations of the CUTs. This approach [8] was the first 
which realized the computation of bound states in higher 
orders [2,3]. The concept of a similarity transformation 
is indispensible for a conceptually clear computation of 
multi-particle effects [4,15]. 

B. Setup 

In Sect. II we analyse global structural aspects of ef- 
fective operators. The basic prerequisite will be Eq. (3). 
Furthermore, we show that the linked cluster property 
holds. Therefore the effective operators which hold in 
the thermodynamic limit can be computed in finite sys- 
tems. 

Sect. Ill is a preparatory section in which the perturba- 
tive CUT for Hamilton operators of certain kind is con- 
structed. Low-dimensional spin models on lattices are 
among the models which can be treated in this way. 

Sect. IV contains a detailed description of how the per- 
turbative CUT can be extended to transform general ob- 
servables. Series expansions in x for the effective observ- 
ables are obtained which allow to compute the experi- 
mentally relevant spectral functions. So the extension 
from Hamiltonians to observables is an important one. 

The article is concluded in Sect. V. 



II. THE STRUCTURE OF EFFECTIVE 
OPERATORS 

In this section, we assume that we are able to construct 
a mapping such that i? c ff fulfills Eq. (3). The eigen-states 
of the particle number operator Q serve as a basis for the 
Hilbert space of the system. If the mapping is realized 
perturbatively, the matrix elements of H e g and O e g are 
polynomials in x. 

A. The effective Hamiltonian 

1. Global structure 

We will show that H c s can be written as 

H cS = H + H 1 +H 2 +H 3 + ... , (5) 
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where H n is an n-particle irreducible operator, i.e. H n 
measures n-particle energies. Moreover, each thermody- 
namic matrix element of any of the components H n can 
be obtained on finite clusters for a given order in x if the 
original Hamiltonian is of finite range. The components 
H n can be defined recursively in ascending order in n. 

Eq. (5) comprises already a route to determine the 
properties of H c h in a sequence of approximate treat- 
ments. The very first step is to know the ground state 
energy which defines H - The second level is to de- 
scribe the dynamics of a single particle (elementary exci- 
tations) correctly which is possible by knowing Hi. The 
third level is reached if H2 is included which contains 
the information on the interaction of two particles. True 
three-particle interactions are contained in H3 and so on. 
From the generic experience in condensed matter theory, 
the three- and more particle terms can very often be ne- 
glected. So the first three terms in Eq. (5) provide the 
systematically controlled starting point of a broad class 
of problems. 

Let us clarify some notation. We define the following 
eigen-states of the particle number operator Q 



at site i. They are bosonic operators. Their definition 
can be tailored to include a hard-core repulsion between 
the particles to account for the common situation that 
at maximum one of the particles may be present at given 
site i. If the particles have additional internal quantum 
numbers, i.e. if there can be different particles at each 
site, the indices i and j are substituted by multi-indices 
i and j. 

As an example let us consider that there are three kinds 
of particles per sites, but that at maximum one of these 
particles can occupy a given site. Then each site cor- 
responds to a four-level system; the particles are hard- 
core bosons. Such a situation arises in antiferromagnetic 
dimerized spin systems where each dimer represents a 
four-level system. The ground state is the unique sin- 
glet while the three particles are given by the three-fold 
degenerate triplet states. In this case we have the multi- 
indices i = (i,a), where i denotes the site and a takes 
for instance the three values of the S z component a £ 
{ — 1,0,1}. In the local basis {\i, s), \i, — 1), \i, 0), \i, 1}}, 
where s denotes the singlet, the local creation operators 
e] „ are the 4 x 4-matrices 



|0) ground state (particle vacuum) 
\i) state with 1 particle on site i 
\i1i2) state with 2 particles on sites ii and 12 



(6) 



i.e. Q\0) = 0|0), Q\i) = and Q\ij) = 2\ij) and so on. 
These states span the global Hilbert space £ of the physi- 
cal system under study. Dealing with (hard-core) bosons 
1*1*2) and \i2i1) are identical states. This indistinguisha- 
bility causes a certain ambiguity. This ambiguity can be 
remedied for instance by assuming that coefficients de- 
pending on several indices I1I2 ■ ■ ■ i n are even under per- 
mutation of any pair of these indices [16]. For simplicity, 
the ground state |0) is assumed to be unique. 

Let 1Z be an arbitrary operator acting on £ and con- 
serving the number of particles [1Z, Q] = 0. By lZ\ n we 
denote the restricted operator acting on £ n C £ spanned 
by all states with exactly n particles. 

Now we define the operators H n 



Ha := Eol 



^2 '■ tjij2;iii2 e j 1 e j 2 e i 2 e ii 



ui2;ji.j2 



H n := £ 



A i e n ' ' ' e j n e i n " ' ' e ir 



(7a) 
(7b) 

(7c) 



(7d) 



where 1 is the identity operator. Note that these opera- 
tors arc defined on the full Hilbert space £. The operators 
e'P are local operators that annihilate (create) particles 



eL 1 = 



-i,0 



e i,i — 













1 























Vo 








0/ 


/° 



















°A 


1 











Vo 






























(! 








:) 










0/ 



(8a) 



(8b) 



(8c) 



It is understood that the action at all other sites but i is 
the identity so that the operators in (8) are defined on the 
whole Hilbert space. The annihilation operators ej jQ are 
given by the hermitian conjugate matrices. All possible 
commutators can easily be computed within the matrix 
representation. Finite matrix elements in the lower right 
3x3 block can be viewed as combined annihilation & 
creation processes: The matrix M a ^ with all elements 
zero except the one at (a, (3) corresponds to the process 
ej a e- q. A finite matrix element in the upper left lxl 
block, i.e. the singlet-singlet channel, can be expressed in 
normal-ordered fashion as I4 — e l a e i a- ^ n this wa y 
the operators (8) and their hermitian conjugate define a 
complete algebra which in turn enables us to classify con- 
tributions of the Hamiltonian according to the number of 
particles affected as done in Eqs. (5) and (7). 

The decomposition (5) is physically very intuitive. Yet 
the next important question is whether and how the op- 
erators H n are unambiguously defined. This issue is ad- 
dressed by noting that H n \ m vanishes for m < n. This 
follows directly from the normal-ordering of the creation 
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and annihilation operators in Eq. (7). Then we can pro- 
ceed iteratively by requiring that H c g applied to n parti- 
cles corresponds to H n + Hi + . . . + H n (n arbitrary but 
fixed). Solving for H n yields the recursions 



#00 


■= #eff|o 




(9a) 


#i|i 


:= #cff 1 


— #o l 


(9b) 


#2 2 


■= #eff 2 


— #02— #12 


(9c) 


#n |n 


:= #cff n 


n-1 
- ^ #i|n . 


(9d) 






i=0 





Assuming that H e g is calculated beforehand one starts 
by evaluating E by means of the first definition. The 
result entirely defines Hq. The restriction Ha\i is then 
used in the second equation to extract the tj-i of Hi and 
so on. Generally, H n is defined on the full many-particle 
Hilbert space, not only for n particles. But it is sufficient 
to know the action of H n on the subspace of n particles 
to determine all its matrix elements in (7). It is the 
essential merit of the notation in second quantization (7) 
that it provides the natural generalization of the action of 
a part of the Hamiltonian on a finite number of particles 
to an arbitrary number of particles. Since Eq. (9d) holds 
for any number of particles and since H n \ m vanishes for 
m < n we obtain Eq. (5), neglecting the precise definition 
of convergence which is beyond the scope of the present 
paper. 

In conventional many-body language, H n stands for 
the n-particle irreducible interaction. The subtractions 
in Eq. (9) ensure that H n contains no reducible con- 
tributions, i.e. contributions which really act only on a 
lower number of particles. It should be emphasized that 
the formalism above does not require that a simple free 
fermionic or bosonic limit exists. It is possible to start 
from any type of elementary particles counted by some 
operator Q. 

Moreover, the formalism presented in this section does 
not depend on how H e g is obtained. It does not matter 
whether a pcrturbativc, a rcnormalizing procedure or a 
rigorously exact method was used to obtain H e g. 



2. Cluster additivity 

Here we focus on formal aspects of a perturbative 
approach generalizing results obtained previously for 0- 
particle properties [17] and for 1-particle properties [18]. 
The feature that the Hamiltonian is of finite range on the 
lattice is exploited. Then the Eqs. (9) can be evaluated 
on finite subsystems (clusters, see below). Still, the ther- 
modynamically relevant matrix elements of the operators 
H n are obtained as we show in the following paragraphs. 

To proceed further definitions are needed. A cluster C 
of the thermodynamic system is a finite subset of sites of 
the system and their linking bonds. By TZ C we denote 



an operator which acts only on the Hilbert space £ c of 
C. If C denotes the sites of the total system which are 
not included in C, the restricted operator TZ C is lifted 
naturally to an operator TZ in the total Hilbert space 
£ = £ c ® £~c by 



TZ := TZ C 



(10) 



Note that it is not possible to define a restricted operator 
TZ C from an arbitrary operator TZ acting on £ since TZ will 
not have the product structure (10) in general. 

Two clusters A and B are said to form a disconnected 
cluster C = AuB iff they do not have any site in common 
A fl B = and there is no bond linking sites from A with 
sites from B. Otherwise the clusters A and B are said to 
constitute together a linked cluster C = A U B. Given a 
disconnected cluster C = AU5 an operator TZ C is called 
cluster additive iff it can be decomposed as 



TZ C = TZ A 



1 b + 1 a 



'TZ' 



(11) 



With these definitions we show that H e g and H n are 
cluster additive. But # e ff| n is n °t! This will turn out to 
be another important reason to introduce the H n . 

The cluster additivity of H^ s is obvious since A and B 
are assumed to be disconnected. So they can be viewed 
as physically independent systems. Hence 



#cff 



H, 



off 



#f ff 



(12) 



Similarly, we deduce from (9) the operators H A and 
which act on £ A and £ B , respectively. Then it is straight- 
forward to verify that the operators 



H? = H A 



1 B + 1 A 



Hi 



(13) 



fulfil the recursion (9) for the operators defined for the 
cluster C. Hence the operators H e g and H n are indeed 
cluster additive. 

It is instructive to see that # e ff|n is not cluster addi- 
tive, contrary to what one might have thought. Let us 
consider the tentative identity 



H 



off I" 



= H 



eff In 



1 B + 1 A 



H 



B I 
off In 



(14) 



This equation cannot be true since on the left hand side 
the number of particles is fixed to n while on the right 
hand side the number of particles to which the identi- 
ties 1 A and I s are applied is not fixed. So no cluster 
additivity is given for the H e g\ n . 

The fact that cluster additivity holds only for partic- 
ular quantities was noted previously for n = 1 [18]. For 
n = 2, the subtraction procedure was first applied in 
the calculations in Ref. [2] (though not given in detail). 
In Refs. [3,4,15,22] the subtractions necessary to obtain 
the irreducible 2-particle interaction were given in more 
detail. The general formalism presented in this article 
shows on the operator level why such subtractions are 
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necessary and where they come from. Thereby, it is pos- 
sible to extend the treatment to the general n-particle 
irreducible interaction. 

The notation in terms of second quantization (7) ren- 
ders the cluster additivity almost trivial. This is so since 
the creation and annihilation operators are defined lo- 
cally for a certain site. It is understood that the other 
sites are not affected. Hence the same symbol ej can be 

used independent of the cluster in which the site i is em- 

t c 

bedded. In particular, one identifies automatically e]' 



ith 



<g> 1 B if i e A and with V 



if i G B. 



Hence cluster additivity is reduced to trivial statements 
of the kind that 



i,j£A 



ijeB 



implies 



i,jeA 



H A <g> 1 B + V 



i,j£B 

>H{ 



(15a) 
(15b) 

(16a) 
(16b) 
(16c) 



In this sense, the notation in second quantization is the 
most natural way to think of cluster additivity. 

Following Gelfand and co-workers [1,17,18] we con- 
clude that the cluster additive quantities possess a cluster 
expansion. Hence all the irreducible matrix elements tj-i 
possess a cluster expansion and can be computed on finite 
clusters. 



3. Computational aspects 

Since H e g conserves the number of particles, i.e. 
Eq. (3), its action is to shift existing particles. Let us 
denote the relevant matrix elements for a linked cluster 
A by 



4i 



-effl 



:= (jlj2|# e ffKl*2) 



(17a) 
(17b) 
(17c) 



where the indices . . . may be multi-indices from now 
on. Put differently, E A is the matrix element of H^\q, 
the a A ^ are the matrix elements of H^\i, the o-^i\i\ii 
those of H A S \ 2 and so on. The number E A is the ground 
state energy of cluster A. The recursive definitions (9) 
imply 



t A - = a A - 

t A ■ =a 

JU2;»ii2 



E A Sji 

I, /_.:/ i 1-2 E Sj 1 i 1 5j 2l2 



(18a) 



-t A 5- -t A S ■ 



t A 8- 



Eo^jii2^j2il 

-tf uil 5 j2i2 (18b) 



a 3lj2j3-Mi2i3 ^0 Al A 2 



(18c) 



where Aq comprises six terms resulting from Hq, A\ com- 
prises 18 terms resulting from Hi and A 2 comprises 36 
terms resulting from H 2 ■ The explicit formulae are given 
in Appendix A. The recipe in deriving the above equa- 
tions is straightforward. For a given n-particle process 
{i m } — ► {jm} {m G {1, . . . ,n}) one has to subtract all 
possible processes which move less than n particles. Since 
the m-particle processes with m < n have been com- 
puted before the procedure is recursive. Note that all 
coefficients must be computed for the same cluster. 

The cluster additivity or, equivalently, the existence of 
a cluster expansion can be exploited to compute the irre- 
ducible matrix elements on finite clusters given that the 
Hamiltonian is of finite range. There are two strategies 
to do so. 

The first strategy is to choose a cluster large enough to 
perform the intended computation without finite-size ef- 
fects. This strategy works particularly well if the dimen- 
sionality of the problem is low. Let us assume for sim- 
plicity that the Hamiltonian links only nearest-neighbour 
sites. Aiming at a given matrix element, for instance 
tf 1 j 2 - 1 i 1 %2i which shall be computed in a given order k, the 
large enough cluster C; contains all possible subcluster 
C s with two properties: (i) they have k or less bonds, (ii) 
they link the concerned sites ji,j2,*i,*2 among them- 
selves [19]. Clearly, C'i depends on the order k. But it 
depends also on the sites j\, j 2 , h, i 2 under study so that 



the notation c[ k \{j\, j 2 ,ii,i 2 }) is appropriate, 
that the order of the sites does not matter. 

If some sites are omitted the constraints for the 
subclusters C s are diminished since less sites must be 
linked. This implies in particular Cj k \{ji, j 2 ,ii,i 2 }) C 
c[ k \{ji,i\\). Hence there can be a cluster A 
which contains Cj k \{ji, j 2 ,ii,i 2 }) but does not contain 
Cj {{ji, h}) so that the hopping matrix element tf ^ is 
not the thermodynamic one, but the interaction t A j 2 -i 1 i 2 
is without finite-size correction. So intermediate steps 
in the calculations (18) can display finite-size effects al- 
though the final result does not. In Refs. [2,3,8,13] we 
followed this strategy. 

The second strategy is to compute for a given order 
k the net contributions of all clusters C with m < k 
bonds which link the sites under study. The advantage 
of this approach is that only smaller clusters need to be 
treated (< k bonds). The price to pay is an overhead in 
determining the net contribution. This requires to deduct 
from the total contribution of C the contributions of all 
subcluster of C with less bonds which link the points 



Note, 
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under study. This must be done in order to avoid double 
counting. More details on this strategy can be found in 
Rcf. [1]. 

For Hamiltonians with relatively simple topology, the 
second strategy is more powerful. For more compli- 
cated Hamiltonians, however, the task to implement the 
overhead without flaw can quickly become impracticable 
while the first strategy can still be used, at least up to a 
certain order of the perturbation. 



B. Effective observables 

An effective Hamiltonian conserving the number of 
particles is useful to determine characteristic energies of 
the considered systems. But it is not sufficient to de- 
termine physical quantities which require more knowl- 
edge than the eigen-energies of the system. In particu- 
lar, we aim at determining dynamic correlations such as 
(O(t)O(0)). Then the mapping of the original Hamil- 
tonian H to the effective Hamiltonian H e s must be ex- 
tended to a mapping of the original observable O to the 
effective observables O e ff- Here we will assume that this 
has been achieved by an appropriate unitary transforma- 
tion, for instance in a continuous fashion as described in 
the Introduction. 



The additional feature in comparison to Eq. (5) is the 
sum over d. Tab. I sketches the structure of the terms 
appearing in the partial observables 0^ n 



d[ 1 n- 
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2 


3 












-3 











eee 


-2 








ee 
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eWeee 
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e T e T ee 


e^eWeee 
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e+ 


e T e T e 


e T e T e T ee 


e T e T e T e T eee 


2 


eV 


e T e T e t c 


e T e T e T e T ee 


e T e T e t e t e T eee 













TABLE I. List of terms appearing in the partial observ- 
ables Od,n which form together the effective observable O e s 
according to Eq. (20). No prefactors or indices are given for 
clarity. 

Let us assume that we computed O c g by some tech- 
nique, for instance by a CUT. Then the partial observ- 
ables can be determined recursively by 







■= Oeff 0- 


-*0+d 




(21a) 


Od.i\i- 




:=Oeff|i- 


-^1+d 


- Od,o\l- 


,i+d (21b) 






:= 0eff|2- 


-^2+d 


- Od,0 2- 


-^2+d - Od,l 2->2+d 



1. Global structure 

The structure of the observables can be described best 
by using the notation of second quantization. Thereby it 
can be denoted clearly how many particles are involved. 
The most important difference compared to the Hamil- 
tonian is that there is no particle conservation. Gener- 
ically an observable creates and annihilates excitations, 
i.e. particles. Hence we define the operators 



Od.n := 



(19) 



E 



w 



J i ■■J n +dlil---i n e j 1 



■]l---3n + d 



The local operators have been described after Eq. (7). 
Again they shall appear normal-ordered, i.e. all creation 
operators are sorted to the left of the annihilation oper- 
ators. The first index d indicates how many particles are 
created (d > 0) or annihilated (d < 0) by application of 
Od.n- The second index n > denotes how many parti- 
cles have to be present before the operator Od, n becomes 
active. The result of Od.n acting on a state with less than 
n particles is zero. 

In analogy to Eq. (5) the effective observables can be 
decomposed into partial observables like 



Coff = Od.r, 



(20) 



n— d>—n 



&d,n \ n- 



>n+d 



O c ff\ r 



i=0 



i \n— m+d 



(21c) 



Here \ n ^ n +d denotes the restriction of an operator to 
act on the n-particle subspace E n (domain) and to yield 
states in the (n+d)-particle subspace £ n +d (co-domain). 
The recursion is set-up in analogy to (9) . It is again used 
that an operator Od, n effectively vanishes if it is applied 
to less than n particles. Barring possible problems to 
define convergence, the validity of the recursion (21) for 
all d and n implies the decomposition (20). 

As for the Hamiltonian the partial observables Od, n 
can be viewed as the n-particle irreducible part of the 
particular observable. The notation in second quantiza- 
tion elegantly resolves the question how the observables 
act on clusters as was explained in the section II A 2. 
Hence the definition (19) ensures cluster additivity and 
there exist cluster expansions for the partial observables. 
So they can be computed on finite clusters. 

If dynamical correlations at zero temperature T = 
shall be described, the observables are applied to the 
ground state |0) which is the particle vacuum [6]. Then 
only the partial observables Od,o with d > matter. 
According to (21a) no corrections are necessary, i.e. the 
structure of the relevant part of the effective observable 
is given by 



Oj ff =° = 0o,O + 01,0 + 02,0 + 03,0 + • • • 



(22) 



G 



This structure has been used so far in a number of investi- 
gations of spectral weights [23,24] and spectral densities 
[13,14,25]. It turned out that it is indeed sufficient to 
consider a restricted number of particles [13,14,24]. But 
the question how many particles are required to describe 
a certain physical quantity sufficiently well depends on 
the considered model, the chosen basis (What do we call 
a particle?) and the quantity under study. 

At finite temperatures a certain number of particles 
will already be present in the system due to thermal fluc- 
tuations. Then the action of the partial observables Od, n 
with n > 1 will come into play as well. This constitutes 
an interesting route to extend the applicability of effec- 
tive models, which were derived in the first place at zero 
temperature, to finite temperatures. 



Pi then the observable O itself links these sites pi. Apart 
from this difference compared to the matrix elements of 
the effective Hamiltonian, we may copy the remaining 
steps from there: 

There are again the two strategies. Either the calcula- 
tion in order k is performed on a cluster C; large enough 
so that all subclusters of k bonds linking the relevant sites 
ji, 32, i,P are comprised in C; [13,14,26,27]. Or one has 
to add the net contributions of all different clusters with 
k or less bonds which link the relevant sites ji,j 2 ,i,p 
[25]. In either way the results for spectral densities can 
be obtained. 



III. TRANSFORMATION OF THE 
HAMILTONIAN 



2. Computational aspects 

The recursive equations for matrix elements which can 
be derived from (21) are very similar to those obtained for 
the Hamiltonian (18). We illustrate this for the matrix 
elements of Oi, n - Let the bare matrix elements on a 
cluster A be 



:= (j\O^\0) 

(hh\o^\i) 



(23a) 
(23b) 



From (21) we obtain the irreducible elements as 



w 



JU2-1 



3 



w h S j2i 



(24a) 
(24b) 



As for the irreducible interactions the strategy is straight- 
forward. One has to subtract from the reducible n- 
particle matrix elements v A the contributions which come 
from the m-particle irreducible matrix elements w A with 
m < n. With this strategy also other irreducible matrix 
elements can be determined in a straightforward manner. 

So far our considerations were general in the sense that 
it did not matter how we achieved the mapping. Next we 
focus on the actual perturbative evaluation of the matrix 
elements on finite clusters. For simplicity, we assume as 
before that the perturbative part of the Hamiltonian links 
only nearest-neighbour sites. Let us consider for instance 
wf 1 j 2 . i - We assume that the observable O is also local, i.e. 
acts on a certain site only, or is a sum of such terms. If the 
observable is a sum of local terms then the transformation 
of each term separately and subsequent summation yields 
the result. So without loss of generality we consider O to 
affect only site p. Then we have to compute the matrix 
elements for clusters linking the four sites ji,j2,i,P- If 
O itself is a product of operators affecting several sites 



So far no particular property of the transformation pro- 
viding the effective operators H e g and O e ff was assumed. 
The only prerequisites were the existence of a counting 
operator Q, which counts the number of elementary exci- 
tations, i.e. particles, and the conservation of this number 
of particles by H e g : [H e g, Q] = 0. 

Here we specify a particular transformation leading to 
[H e g,Q] = 0. This section is a very brief summary of 
Ref. [8] which is necessary to fix the ideas and the nota- 
tion for the subsequent section dealing with the transfor- 
mation yielding the effective observables. 

For simplicity we restrict the considered systems in 
the following way: The problem can be formulated as 
perturbation problem as in Eq. (2) with the properties 

(A) The unperturbed part U has an equidistant spec- 
trum bounded from below. The difference between 
two successive levels is the energy of a particle, i.e. 
Q = U. 

(B) There is a number N 3 N > such that the perturb- 
ing part V can be split as V = X^=-iv ^™ where T n 
increments (or decrements, if n < 0) the number of 
particles by n: [Q, T n ] = nT n . 

Condition (A) allows to introduce the particularly sim- 
ple and intuitive choice Q = U. Note that the restric- 
tions of (A) are not too serious in practice since very 
often the deviations from an equidistant spectrum can 
be put into the perturbation V. Conditions (A) and (B) 
together imply that the starting Hamiltonian H has a 
block-band-diagonal structure as depicted in Fig. 1. The 
perturbation V connects states of different particle num- 
bers only if the difference is a finite number < N. Note 
that very many problems in physics display this property, 
for a discussion of interacting fermions see Ref. [5,6]. So 
far, most applications consider N = 1 [3,28] and N — 2 
[2,8,10,13,14,21,24,26,27,29], but calculations for higher 
N are also possible [30]. 

We solve the flow equation (4) for the Hamiltonian (2) 
obeying the conditions (A) and (B) pcrturbatively, that 
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means up to a certain order in the expansion parameter 
x. The ansatz used is 

oo 

H(x;£) = U + ^x k F (^ni)T{m), 

k=l |m|=fc (25) 

with unknown real functions F(£; m) for which the flow 
equation (4) yields non-linear recursive differential equa- 
tions [8]. The notation comprises 



m = (mi , 77i2 , fi3 , • • . , TO/t ) with 
mi e {0,±1,±2,... ,±N} 
\ m \ = k 
T ( tti ) — T mi T m2 T m3 • • • T 7nk 



M(m) = 



(26a) 
(26b) 
(26c) 
(26d) 

(26e) 



The second sum in ansatz (25) runs over all indices m 
of length |m| = k. Thereby, H(x; £) includes all possi- 
ble virtual excitation processes T(m) in a given order x k 
multiplied by the weight F(£; m). 

The optimum choice for the infinitesimal generator r\ 
of the unitary transformation reads 

oo 

r)( X ; £) = Y x k Y sgn (M(mj)F(£]m)T(m). 

k=l \m\=k (27) 

In the eigen-basis {\n}} of Q, i.e. Q\n) — n\n), the matrix 
elements of the generator r\ read 



r) ld (x;£) = sgn(<3i - Q j )H iij (x;£) 



(28) 



with the convention sgn(O) = 0. This choice keeps the 
the flowing Hamiltonian block-band diagonal also at in- 
termediate values of £ [8,11]. For £ — > oo the generator 
(28) eliminates all parts of H(x;£) changing the number 
of particles so that [H c g , Q] = with H c r :— H{£ = oo). 

For the functions F(£;m) a set of coupled differential 
equations is determined by inserting Eqs. (25) and (27) 
in the flow equation (4) and comparing coefficients. The 
differential equations are recursive [8]. The functions F 
of order k + 1, i.e. F(£;m) with \m\ = k + 1, are deter- 
mined by the functions F of order k. The initial con- 
ditions are F(0;m) = 1 for |m| = 1 and F(0;m) = 
for \m\ > 1. The functions are sums of monomials with 
structure (p/q)£ l exp(— 2/j£), where p, q, i, (p > 0) are in- 
tegers. This allows to implement a computer-aided iter- 
ative algorithm for the computation of the functions F 
[8]. 

The following symmetry relations hold 

F(£;m) = F(£; (-m k , . . . ,-m x )) (29a) 
F{£-m) = F(£- {-m x , . . . , -TO fe ))(-i)'- l+1 . (29b) 

Relation (29a) reflects the hermitecity of the Hamilto- 
nian. The block-band diagonality for all £ implies 



In the limit £ — > oo the coefficients C(m) := F(oo;m) 
are obtained. They are available in paper form [8,28] and 
electronically [31]. The effective Hamiltonian is given by 
the general form 



H eS (x) = U + 



E* fc 

fc=i 



E C (m)T(m), 



2i) =0 



(31) 



where M (m) = reflects the conservation of the num- 
ber of particles. The action of i? e ff can be viewed as a 
weighted sum of particle-number conserving virtual exci- 
tation processes each of which is encoded in a monomial 
T(mj. We want to emphasize that the effective Hamil- 
tonian H e g with known coefficients C(m) can be used 
straightforwardly in all perturbative problems that meet 
conditions (A) and (B). 





Hoi 








H n 








H21 




H23 






H32 


H 33 



FIG. 1. Block-band diagonal Hamilton matrix for N = 1 in 
the eigen-basis {|n)} of the operator Q which counts the num- 
ber of particles. The unperturbed Hamiltonian H(x — 0) =U 
and the effective Hamiltonian H c ff have matrix elements in the 
dark areas only: [H c s,Q] = 0. For non-degenerate ground 
state -ffoo is a 1 x 1 matrix. The dimension of H„ n grows 
roughly like L" with system size L. The perturbation V can 
lead to overlap matrices indicated as light boxes. The empty 
boxes contain vanishing matrix elements only. 



IV. TRANSFORMATION OF OBSERVABLES 

To calculate physical quantities which do not only de- 
pend on the eigen-energies the relevant observables must 
also be known. The conceptual simplicity of unitary 
transformations implies that the observables must be 
subject to the same unitary transformation as the Hamil- 
tonian. In this section we describe how the perturbative 
CUT method can be extended to serve this purpose. 

Consider the observable O. It is mapped according to 
the flow equation 



F(£;m) = for \M{m)\> N 



(30) 



dO(x;t) 
d£ 



[r){x;£),0{x;£)] 



(32) 
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where the same generator r](x; £), given in Eq. (27), as in 
Eq. (4) is to be used to generate the transformation. In 
analogy to Eq. (25) we employ the ansatz 

oo fe+1 

k=0 J = l \m\=k (33) 

where the G(£;m;i) are real- valued functions for which 
the flow equation (32) yields recursive differential equa- 
tions. The operator products 0(m;i) are given by 



0(m-i) :=T mi ---T mi _ 1 OT„ 



(34) 



where we use the notation of the Eqs. (26). The integer 
i denotes the position in Q(m, i) at which the operator 
O is inserted in the sequence of the T m . The starting 
condition is 0(x; 0) = 0{x) and the final result is found 
at £ = oo: O cS (x) := 0(x; oo). 

Inserting the ansatz (34) for 0(x;£) and the generator 
r](x; £) from (27) into the flow equation (32) yields 



fc+i 



c) 



k=0 \m\=k i=l 
oo oo 



fe 2 +l 



E E^ 1+fc2 E E^^)G(^m'V) 



\rn,'\ = k 1 1=1 
\m"\ = k 2 



fci=l fc 2 =0 

xsgn(M(m')) [T(m'), 0(m"; i) 



(35) 



The functions F(t; to) are known from the calculations 
described in the previous section III pertaining to the 
transformation of the Hamiltonian. The sums denoted 
by expressions of the type |m| — k run over all multi- 
indices m of length k. 

Comparing coefficients in Eq. (35) yields a set of re- 
cursive differential equations for the functions G(£;m,i). 
To ease the comparison of coefficients we split a specific 
to with k fixed in two parts as defined by i 



m = (mi, m r ) , 



(36) 



with |to ; | =i — l and \m r \ = k — i + 1 such that the 
splitting reflects the structure of 0(m;i) in (34). Then 
the exlicit recursions can be denoted by 

^G(£;to;z) = 

sgn(M(TO a ))F(£;TO a )G(^; (m b ,m r );i- |toJ) 

- s MM(m b ))F(£; mb )G(£;( mi ,m a );z) . (37) 

nh r = {m_ a ,m_ b ) 

The recursive nature of these equations becomes appar- 
ent by observing that the summations m ; = (m a , m b ) 
and TO r = {m al m b ) are performed over all non-trivial 
breakups of to ; and to,.. For instance, the restriction 



TO ; = (mi, m2, ■ ■ ■ , m i-i)==(zB;a: Hh>) w hh m a means, 
that one has to sum over the breakups 



TO a = (mi) 

TO a = (toi,TO 2 ) 



and m b = (m 2 , . . . ,mj_i) 
and m b = (m 3 , . . . ,mj_i) 



rn a = (toi,to 2 , . . . ,to;_i) and m b = () . 



(38) 



This implies that the G(£; to; i) appearing on the right 
side of Eq. (37) are of order k — 1 or less. Once they are 
known the function on the left hand side of order k can be 
computed. By iteration, all functions can be determined. 
The initial conditions follow from 0(x; £ = 0) = O and 
read 



G(0;to;1) = 1 for |m| =0 
G(0;m;i) = for \m\ > . 



(39a) 
(39b) 



By iteration of (37), all functions can be determined. 

We briefly discuss two examples to illustrate how the 
Eqs. (37) work. Let us assume N = 2. All zero order 
functions G(£;(),l) are equal to 1. Since there is no 
breakup of (), as would be required by the sums on the 
right hand side of Eqs. (37), the right hand sides vanish 
identically, whence G(£; (); 1) = 1 for all values of £. 

The first order function G(£; (1); 2) is given by 

^G(£; QJ; 2) = sgn [M((l))] F(£; (1)) • G(£; (); 1) 



= e 



(40) 



where F(£; (1)) = e £ is taken from Eq. (15) in Ref. [8]. 
With the initial condition G(0; (1); 2) = from (39) the 
differential equation (40) yields 



G(4(l);2) = l-e- 



^oo 



(41) 



As a second example we consider a second order func- 
tion where we can use the above result 



d_ 

d£ 



G(^;(^2 j _l);3) = 

mi 

sgn[M((-2,l))]F(^(-2,l))G(£; (), 1) 
+sgn [M((-2))] F(£- (-2))G(£; (1), 2) (42a) 
= - (e- 3£ - e- e ) ■ 1 - e- 2e ■ (l - e~ l ) (42b) 
= e~ l - e~ n . (42c) 



Again the functions F are taken from Eq. (15) in Ref. [8]. 
Integrating the result (42c) using the initial condition 
(39) leads to 



G(£;(-2,l);3) 



+ h~ 2e 



l— »oo 



(43) 
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This kind of calculation carries forward to higher orders. 
The functions G - like the functions F - are sums of sim- 
ple monomials {p/q)t exp(— 2/j£), where p,q,i,(/i > 0) 
are integers. Thus the integrations are always straight- 
forward 



/' 

Jo 



d£'. 



1 ei+l 



/' 

Jo 



dl'te 



2\i 



(2 M )< 



3=0 



(2 M )' 



(44a) 
(44b) 



and can easily be implemented in a computer-algebraic 
programme. The remaining implementation follows very 
much the same line as described previously for the func- 
tions F [8]. 

In analogy to Eqs. (29) for F two symmetry relations 
hold for G. With m — (mi, . . . , mk) they read 

G(£; m; i) = G(£; {-m k , ... , -mi); k-i + 2) (45a) 
G(£;m;i)=G(^;(-mi,... , -m k ); (45b) 

as can be shown by induction. The first symmetry (45a) 
holds if O is hermitian. Unfortunately, there is no equiv- 
alence to Eq. (30) so that a possible initial block-band 
structure in 0(x; 0) is generically lost in the course of 
the transformation, i.e. for £ > 0. 

In the limit £ — > oo the coefficients C(m:i) := 
G(oo;m;i) E Q arc obtained as rational numbers. So 
we retrieve finally 

oo fe+1 

O cS (x) =J2 xk J2 C(m;i)0(m;i) 

k=0 i=l \m\=k (46) 

similar to Eq. (31). We will make the coefficients C(m; i) 
available electronically [31]. Note that O c g is not a 
particle-conserving quantity as is obvious from the fact 
that the sum over |m| is not restricted to M(m) = 0. In 
order to see the net effect of O e s(x) on the number of 
particles explicitly it is helpful to split the bare opera- 
tor accordingly O — X^=-JV' where T' n increments 
(or decrements, if n < 0) the number of particles by n: 

[Q' T n] = nT n- 

The difference between the bare initial observable O 
and the representation (46) must be viewed as vertex cor- 
rection which comes into play since the bare initial excita- 
tions are not the true eigen-excitations of the interacting 
system. We like to stress that the formalism presented 
introduces the notions of n-particle irreducibility, vertex 
correction and so on without starting from the limit of 
non-interacting conventional particles such as bosons or 
fcrmions. 



V. CONCLUSIONS 



A. Summary 



In this article we have presented an approach to calcu- 
late energies and observables for quantum multi-particle 
systems defined on lattices. The article has two main 
parts. In the first part (Sect. II), we assumed the exis- 
tence of a mapping of the original problem to an effec- 
tive one in which the number of elementary excitations, 
the so-called (quasi-)particles, is conserved. The general 
structure of the effective Hamiltonians and the observ- 
ables is analysed. We found that a classification of the 
various contributions in terms of the number of parti- 
cles concerned is most advantageous. To this end we 
introduced a notation in second quantization which does 
not, however, require non-interacting fermions or bosons. 
Generically, hard-core bosons are involved. 

We found the formulation in second quantization par- 
ticularly intuitive. It provides in a natural way the ir- 
reducible quantities on the operator level which display 
cluster additivity. We like to emphasize that the defi- 
nition of irreducible operators is not a trivial task if a 
strong-coupling situation is considered as was done in 
the present paper. No limit of non-interacting bosons or 
fermions is assumed. Since the definition of irreducible 
operators is completely general it allows to compute the 
n-particle contribution for arbitrary n. For instance the 
formulae for the 3-particle interactions are given for the 
first time in the literature. 

The irreducible interactions and vertex corrections pos- 
sess a cluster expansion so that they can be computed on 
finite clusters provided that the Hamiltonian is of finite 
range. This property is the basis for the real-space treat- 
ment of many spin systems. 

In the second part (Sects. Ill and IV), we described an 
actual mapping which provides effective operators. The 
mapping is based on continuous unitary transformations. 
In this paper we constructed the mapping perturbatively 
(see 'Outlook'). In Sect. Ill the treatment of the Hamilto- 
nian is given. The computation of the effective Hamilto- 
nian requires the solution of a set of recursive non-linear 
differential equations. For the perturbative set-up under 
study these equations can be solved in full generality, i.e. 
no particular details of the model must be known. 

In Sect. IV we have given the calculational steps to 
compute effective observables. Again recursive differen- 
tial equations have to be solved. But they are linear since 
the transformation of the Hamiltonian is known. For the 
perturbative set-up under study also the equations for 
the observables can be solved in full generality, i.e. no 
particular details of the model must be known. 

The above approach has been used to compute spec- 
tral functions, i.e. dynamical correlations, in a number of 
models [2,3,8,10,13,14,21,24,26-30]. These results may 
serve as examples for the utility of the approach pre- 
sented. 
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B. Outlook 

We like to point out two important consequences of the 
formulation of the effective operators in second quanti- 
zation. Both implications are based on the observation 
that the irreducible operators are defined on the whole 
Hilbert space, i.e. not only for a small number of parti- 
cles. The matrix elements of the n-particle irreducible 
operators can be computed considering only n-particles. 
But the resulting operators hold for arbitrary number of 
particles. 

a. Consequence 1: The effective Hamiltonian is valid 
at finite temperatures. Hence it is possible to extend the 
results obtained in the first place at zero temperature to 
finite temperatures. The technical difficulty arising is to 
treat the interactions properly, in particular the hard- 
core constraint. But the description in terms of effec- 
tive particles helps to tackle this situation. Let us recall 
that at zero temperature no excitation, i.e. no particle, 
is present. At low temperatures only a small density of 
particles will be in the system. So it is well justified to 
use a ladder approximation. This approximation is also 
suited to deal with the hard-core constraint (Bruckner 
approach) [32,33]. Note that the problems linked to the 
existence of anomalous Green functions [33] do not occur 
if the particle-conserving, effective Hamiltonian is used. 
Thus the Bruckner approach for the effective Hamilto- 
nian after a suitable mapping [4,8] is well justified and 
represents a very promising route to treat finite temper- 
atures. 

b. Consequence 2: So far the mapping to an effective 
model has been constructed perturbatively. That means 
that all operators, the Hamiltonian H, the generator rj 
and the observables O, are given in a series of some small 
parameter x. In actual applications these series are suit- 
ably extrapolated. But a certain caveat persists if the 
starting point is a local Hamiltonian. Then a calculation 
up to a certain order describes processes of a certain finite 
range only. This restriction can be partly overcome by 
extrapolating in momentum space, e.g. for dispersion re- 
lations uj(k). But it is difficult to extrapolate the matrix 
elements of the 2-particle irreducible interaction because 
it is not diagonal in all momenta. 

This problem can be overcome by performing the con- 
tinuous unitary transformation directly on the level of 
the n-particle irreducible operators. An ansatz for the ef- 
fective Hamiltonian is chosen comprising for instance all 
possible irreducible n-particle terms and similar terms 
creating and annihilating particles. This ansatz is in- 
serted in the flow equation (4). Comparison of the co- 
efficients tj 1 _..;i 1 _.. and ditj 1 ...-i 1 ... in front of the terms 
. . . e it . . . yields coupled non-linear differential equa- 
tions. These differential equations represent renormaliza- 
tion equations for the problem under study. We call this 
type of transformation a self-similar one since the kind 
of terms retained stays the same. Again, the formulation 
in second quantization allows a significant generalization. 



We like to stress again that the approach presented does 
not require a weak-coupling limit. For illustration, how- 
ever, the reader is referred to the weak-coupling examples 
in Refs. [5-7]. 

c. Concluding Remark In this paper we discussed the 
general structure of effective operators and a perturba- 
tive unitary transformation to derive them. No concrete 
application is presented since such applications can be 
found elsewhere. The two continuative points above show 
along which lines the general structure can be exploited 
to extend the applicability beyond zero temperature re- 
sults and beyond finite range processes. Work along these 
lines is partly under way, but deserves definitely further 
attention. 
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APPENDIX A: THREE-PARTICLE 
IRREDUCIBLE INTERACTION 

Here we complete the formulae for the irreducible 3- 
particle interaction which was given in Eq. (18). The 
corrections Aq, A\ and A2 result from Hq, Hi and H2, 
respectively, as given in (9). They read 



Aq — E [Sj 1 i 1 5j 2 j 3 -i 2 i 3 + Sj 1 i 2 Sj 2 j 3 -i 1 i 3 + Sj 1 i 3 Sj 2 j 3 jtiij] 



(Ala) 







Yt A ■ 


"32j3',ili3 


Vt A ■ 


$ 32 33 


;*i«2 


t A . 

J2;«i 




vt A ■ 

r L 32-,l2 


3ji33;iii3 ' 


Vt A ■ 

r 1 32^3 


^3133 


\ili2 


t A ■ 

33\n 




^t A ■ 


^3i32;iii3 


Vt A 

r l 33V3 


^3132 


;iii2 



M — $jliltf 2 j 3 ;i 2 i 3 + ^jl»2*jtj 3 ;iii3 + &3H3 t f 2 j 3 :i 1 i 2 + 



3233 

ZA 



r 7 A 1 c 7 A 1 c 7 A 

°32il l j 1 j 3 ;i 2 i3 hi2 l j 1 j 3 ;i 1 i 3 32i3 l j 1 33;i 1 i 2 



3233Y 



^'3*i*jij 2 ;»2*3 + ^'3»2*jij 2 ;iii 3 + ^33i3 t j 1 j 2 ;ili2 (^ 1C ) 



(A2a) 



where we used the shorthands 

$3l32\ili2 ■— 3jlil^32i2 + ^jit2^j'2»i 
/ 1 . . • / 1 . . + f A . ■ - + t A ■ ■ ■ + t A ■ ■ 

3l32;ill2 ' 3132--1112 1 3l32;l2ll 1 J2ji;»l»2 1 J2ji;*2*l" 

(A2b) 

While the actual formulae are lengthy the underlying 
principle is straightforward (see main text). Note that 
in concrete realizations it is often advantageous to de- 
note only one representative of the states which do not 
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change on interchange of particles = \ij)). Fur- 

thermore, certain problems allow to exploit higher par- 
ticular symmetries like spin rotation symmetry. Then 
additional permutation symmetries among the various 
quantum numbers constituting the multi-index can be 
exploited leading to the appearance of exchange-parity 
factors. 
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